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Abstract 

For ergodic adiabatic quantum systems, we study the evolution 
of energy distribution as the system evolves in time. Starting from 
the von Neumann equation for the density operator, we obtain the 
quantum analogue of the Smoluchowski equation on coarse-graining 
over the energy spectrum. This result brings out the precise notion of 
quantum diffusion. 
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There are many physical situations where a separation of time scales oc- 
curs, the adiabatic approximation to an evolution is one of the expressions 
of such a separation JTJ. Classically, an important problem in this context 
is to describe the evolution in phase space of an ensemble of systems under 
an ergodic adiabatic Hamiltonian . An ergodic adiabatic Hamiltonian can 
be defined as one which describes ergodic (usually specialized to chaotic) 
dynamics under a slowly time-dependent Hamiltonian. In this paper, we 
study the quantal version of this treatment and obtain results which, in an 
appropriate limit, reduce to results in [||. The reasons to present quantum 
analogue of the classical case study are many-folded. Firstly, classical physics 
being a special case of quantum physics, it is interesting to know what is the 
equation that one gets so that the limiting situation is clearly interpreted. 
Secondly, multiple time-scale analysis is usually employed on classical equa- 
tions and it is, therefore, interesting to expose how this method may be 
shown to be fruitful for quantal equations also ||. Thirdly, a clear analysis 
based on time-scale separation of the von Neumann equation paves way for 
the next-order complication that arises in combining the adiabatic and semi- 
classical limits [|J . The motivation behind this combination of singular limits 
is provided by a work where first-order velocity-dependent corrections to the 
lowest adiabatic approximation for the reaction force on the slow system 
are studied [H. In the classical setting, one recovers geometric magnetism 
and deterministic friction as the reaction forces whereas in the half-classical 
setting (fast system treated quantum mechanically) there is only geometric 
magnetism, no friction. Deterministic frictional force was found in |J and is 
non-zero when the fast motion is classical and chaotic. 

Very interesting examples can be cited from different fields of physics and 
chemistry that resonate with the abovementioned problem. One of the well- 
studied problem is when there are many particles moving in a time-dependent 
shape of the box 0. This is an idealization of nuclear fission and fusion. It 
has been shown in a numerical experiment that the transition from ordered 
to chaotic nucleonic motions is accompanied by a transition in collective 
properties of nuclei from those of elastic solid to visco-elastic to viscous fluid 
||. Quantum mechanical origin of dissipation in finite Fermi systems has 
been recently explained |J. It has been found that the geometric phase, 7 
acquired by a single-particle wavefunction in adiabatic perturbation is related 
to the absorptive part of the frequency-dependent response function, x" °f 
the finite bulk, hence dissipation - we call them 7 — \ relations. Mesoscopic 
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systems can also show this kind of behaviour in their conductivity properties. 

We now describe the basic problem. Consider a quantum system evolving 
adiabatically in time. Since the evolution is adiabatic, we have energy levels 
at every instant of time and an instantaneous basis. Let us now assume that 
the energy levels do not cross and the reason for this be left unspecified. 
In particular, this may arise in a random matrix hypothesis for the system. 
During the evolution, although the levels do not cross, they may come arbi- 
trarily close to each other. Very small spacing between levels then leads to an 



increased probability of non-adiabatic Landau-Zener transitions which 



eventually change the energy distribution of the system. The basic idea be- 



hind using the Landau-Zener transitions goes back [TT[ in the literature of 
nuclear physics where it is used to explain damping of collective modes. 

The classical version of the problem treated here has a distinguished his- 
tory. Starting from the earlier works of Ehrenfest, in classical context, it was 



shown first by Lenard |12[ that if Hamiltonian, H =constant for t < and 
t > ii, then the values of the reduced action § pdq for t < and t > ti 
differ from each other by 0(e m ), however large m may be. Generality of 
the adiabatic invariant led Ott |13[ to show that the error in ergodic adia- 
batic invariant in the classical version of our problem is diffusive, however 
the equation differed from the Smoluchowski equation. The problem was 
re-examined by Jarzynski pj where the Smoluchowski equation was restored 



(although the authors of |2| call it by the name of the Fokker-Planck 



equation). The quantum version is the subject here. To start with, due to 
important differences between classical and quantum mechanics, and with 



lack of a precise notion of chaos in quantum systems [|14[] , it is not clear what 
will the final equation for energy distribution be. We now show that the 
equation is different from the classical case. 

To begin the analysis, let us consider the Hamiltonian, 

H(t) =H + etV (1) 

where H , V are linear operators. We assume that the evolution in time is 
adiabatic which corresponds to the smallness of e. As a concrete realization 
where avoided level crossings may occur, one may assume that the linear 
operators H an d V belong to some random matrix ensembles invariant un- 
der a canonical group where it may be mentioned that a detailed study on 
time correlation functions already exists ||15|| . Thence, at any instant, sys- 



tem admits an eigenvalue spectrum given by the eigenvalue problem for the 
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"frozen" Hamiltonian, 

H{et)\n{et)) = E n {et)\n{et)). (2) 

The probability of the system residing in the state \n(et)) at time t is given 
in terms of the density operator, p, 

p n (t) = (n(et)\p\n(et)), (3) 

with 

E Pn(t) = lVt. (4) 
n 

On the other hand, the probability of finding the system in the range AE 
about energy E, is defined through 

p(E)dE = E p n (t) (5) 

E<E n <E+AE 

With ,00 

/ dEp(E) = 1. (6) 

J— oo 

To relate p n and p(E), let us define the density of states, 
so that 

p(E)AE = P(E)^- = p(N)AN (8) 

implying thereby 

p(E) = E(E)p(N). (9) 

Eq. (9) relates the probability of finding N levels in the energy range [E, E + 
AE] and the probability of finding the system in the state corresponding to 
energy E. 

From time t = 0,the levels evolve in time, the resulting density operator, 
p satisfies 

ih^ = [H,p]. (10) 
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Starting from (10), our objective is to derive an equation for the energy 
distribution, 

rE 



dE'tr{6(E' - H)p}. 



(11) 



Since e in (1) is small (adiabaticity) parameter, we have two time-scales - 
t ("fast" scale) and et ("slow" scale). To incorporate these scales in the 
problem, we employ the multiple time-scale method for treating the partial 
differential equation (10). Accordingly, denoting the set of instantaneous 
states by {\n(et))}, we can write an expansion for the density operator, 

p{{\n{et))},t) = p {{\n(et))},et) + ep 1 {{\n{et))},t,et) + ... (12) 

with the initial conditions, 

Po = Poo(H(et)), pi = p2 = Ps = ■■■ = 0. (13) 

Substituting (12) in (10), we get a system of equations separated by different 
orders of e : 



[po,H(et)] = 0, 



-ih 



dps 



j-1 



det 



1,2,... 



(14) 
(15) 



If there are no other constants of the motion than H(et) on the fast scale, or 
under the Thomas- Fermi approximation, by (14), 



p ({\n(et))},et) = p' (H(et),et) 



(16) 



where the arbitrariness of p' is removed by insisting that p remains valid for 
times 0(e _1 ) by removing secularities in (15) with j = 1. To realise this, we 



operate on the j = 1 equation by an arbitrary operator- valued function [16 



g (H) and perform the trace of the resulting equation over the frozen basis, 



n 



dpc 



d{et) 



ii 



(17) 



For po to be valid for times for 0(e x ), the right hand side (RHS) of (17) 
should be set to zero, which leads to 



5 



dp' (H) OH dp' (H) 
dE d(et) d(et) 



n ) = 0. 



Let us now define 

E(£,et) := "£{n\6(E - H)\n) 

n 

d 



dE 



£<n|0(£ - H)\n) := 



ggCg, et) 

as 



also, we define an energy average over the "frozen" Hamiltonian by 



which implies 



Now, (18) becomes 



-J2(n\5(E-H)...\n) := (...) £ , et , 



£>|...|n) = I dEX(...) E , et . 



dp'o I OH 



df>' 



0. 



dE \d{et)/ Eet d{et) 
Calling 

and using (19), we obtain the identity (see Appendix A), 

<9E d 



d(et) dE 



(£u) = 0, 



(18) 



(19) 



(20) 



(21) 



(22) 



(23) 



(24) 



which is a new derivation of quantum adiabatic theorem. Therefore, (22) 
reduces to 

° <Po£) + 4(^Po) = 0. (25) 



d(et) 



dE 



Then, for p , we have 

dpo 
d(et) 



dE' 



dH 



(26) 
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With this equation and the initial condition (13), we have completely deter- 
mined po- 

The notation d x stands for a partial derivative with respect to x. 

We now proceed to determine p\. The formal solution of (15) with j = 1 

is 

Px({\n)},t,et) = pu + pih 

1 dt'^h({\N)},et) + fc(H{et),et), (27) 



o 



d(et) 



where \N) = \N)(\n),t,t',et) is a state reached at time t' evolving back- 
ward the state |n) at a time, t, under H(et). To determine p[, we remove 
secularities at 0(e 2 ) by a similar procedure as above, resulting in 



0_ 

Of 



^2(n\gp 2 \n} = - Wng 



dpi 



•diet) 



ii 



(28) 



or, more explicitly, 

/ dE J t £H*(£ - H)gp 2 \n) = - J dE^(n5(E - H)g 
dEY,(n\S(E - H)gd et f dt'd et p ({\N)},et)\n) 
J dEj2(n\S(E - H)gdMn) :=T\+T 2 . 



dpi 



d(et) 



n 



(29) 



Ti and T 2 are the abbreviations of the two terms above that line in (29). 
Using (26), Ti can be written as 



Ti = J dEj2(n\gS(E - H)d et jf* dt'd E p' (d £t H - u) \n). 



(30) 



We now employ the notion of distributional or weak derivative of distributions 
(denoted by S'(x) in the case od delta distributions). Employing the following 
property of the Dirac delta distributions [17|], viz., 



5(E - H)d et § = 8'(E - H)d et H<!>, 



(31) 
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(30) becomes 

T x = [ dEj2(n\gS'(E - H)d et H f dt'd E p' (d et H - 
J Jo 

= -J dEY,(n\g5{E - H)d 2 E p' d tt H J* dt'(d et H 
- J dEj2(n\gS(E - H)d E p' d et H £ dt'd E (d et H 



u n) 



u)\n) 



u)\n) 



(32) 
(33) 



where use has been made of 

5'(E-H)d E p' = -5(E-H)d 2 E p' . 

This complicated set of terms can be simplified somewhat. To do so, we 
rewrite 7\ in a way that will help us in introducing two-time correlation 
functions later. So, 

T x = -JdE^(n\gS(E-H)dlp' iys(d €t H({\n)})-u)(d €t H 

- J dEY,(n\g5(E - H)d E p' jy s (d et H({\n)})-u)d E (d et H({\N)})-u)\n) 

- J dEj2(n\g5(E - H)d E p' J° ds ud E (d et H - u)\n) 

(34) 

The average two-time correlation function can now be introduced : 
C et (s,E) = ({d et H({\n)},et)-u}{d et H({\N)},et)-u}) E>et 
= W ( n \S(E - H){d et H({\n)}, et) - u}{d et H({\N)}, et) - u}\n)zh) 

Before getting back to (34), we note some simple relations. First of all, we 
can write : 



dsC(s) 



/_/4 S ( n ^ E - H){d f _ t H({\n}}, et) - u}{d £t H({\N)}, et) - u}\ 
f dsJ2(4(E - H){d et H({\n)}) - u}{d et H({\N)}) - u}W) (36) 



l <9£ r° 
E9£ 



Secondly, 



A( s />c W ) = o. 

With (36) and (37), we can now write a relation, 



(37) 



(38) 



which we shall use shortly. 

To simplify T% and T 2 , we have to employ some further averaging pro- 
cedure. We call this a coarse-graining in which we replace a function of an 
eigenvalue, E n by some average quantity such that the explicit dependence 
on the label n disappears. One of the ways it may be done is by an integra- 
tion over the average density of states. The essential point about averaging 
is that the spectrum "seen" by the system is a continuous one. 

After some tedious manipulations, repeated usage of the properties of 
distributions |I7|, and effecting coarse-graining, we arrive at 



Ti = -2 J dEg{E)d 2 E p' n Jl j\sC et (s,E) 



dEg(E)d E p' nn Ed E J^dsC et (s,E) 



o 



dEg(E)d E p' nn {n\S(E - H)d E (d et H({\n)})) - u) 
ds{d tt H{{\N)}))-u)\n) 



dEg(E)d E p' nn u(n\6(E~H) 



0, 



ds(d et H({\N)}))-u) 



\n), 



(39) 



where we have employed ad hoc coarse-graining and replaced 

g(E n ) by g(E), and 



dE 2 



n 



by 



d 2 p' 
dE 2 ' 



(40) 



The last two terms are of the same order and opposite sign, so they will 
simply compensate for each other. 
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Because, for times of 0(j), 

J _dsC(s) = -j_ go dsC(s):=-G 2 , 

we can write for T\ : 



' dE dE' 



(41) 



(42) 



Notice that T 2 has the same form as the expression involving p' in (18), 
manipulations are identical. Finally, the condition that removes secularities 
to C(e 2 ) is 



2 dE dE 



(43) 



In terms of occupation probabilities, po and pi (Cf. Eq. (3)), we have 

(44) 



_( PO E) + _(„ PO E)=0, 



(pi E) + A (upiS) _4( EG2 ^ 



^(et)^ 1 ^ ' dE^'"' dE 
The energy distribution, defined as 

follows the following equation, 

d 



dr) d 2 d 

m = ~ e dE {ur]) + e dE 



G ^dE ' - 



1 d Po dG 2 
2 E dE^E = °" (45) 



(46) 



e 2 <9G 2 g 
+ 2" <9£ dE 



(47) 



which is the final result. This equation is different insofar as there is an 
extra term as compared to the Smoluchowski equation. Thus, the diffusion 
in quantum systems has to be qualitatively and quantitatively different as 
the diffusion coefficient will be different from the one we have in the Smolu- 
chowski equation. 

It is clear that the difference between (47) and the Smoluchowski equa- 
tion is the derivative of tegrated time correlation function. This is, indeed, 
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reminiscent of the relations between friction and diffusion coefficients in weak- 



turbulence plasma theory fL8| . This brings us to the premise on which we 
began, the time scales. 

First of all, the time scale associated with the decay of correlation func- 
tion, 

t c := [C(O)]- 1 f +CO C(s)ds. (48) 

If the quantum system considered is modelled as a random matrix of di- 
mension N [|]1| with large N, we know that correlation function will decay 
very rapidly. Thus, t c can be very small if the quantum systems possess the 
following properties : (a) the number of eigenvalues is very large, and the 
energy spectrum is complex, and, (b) the corresponding classical system is 
chaotic. Chaos in the underlying classical system plays a fundamental role in 
the decay of correlation functions. It was recently shown |TD[ that the quan- 
tum time-depende correlations in a Fermionic system are dominated by the 
classical correlation function. The decay of the correlation function is shown 
in this work to be governed by the eigenvalues of the Liouvillian operator. 
Thus, t c is related to the Liapunov exponents and other detailed features of 
chaos. In classical ergodic adiabatic systems, the time t (fast scale) is much 
larger than t c , thus the third term of (47) is zero. However, in quantal sys- 
tems, we have the quantum mechanical scale, t q = h/S (S being the mean 
level spacing) which is why the third term at 0(e 2 ) is explicitly present. If 
t q <C t c <C t, the quantum effects will dominate, and all the terms in (47) 
will be important. If t c t q t, then the system will behave classically 
initially and eventually, quantum phenomena will become important; so ini- 
tial evolution will be Smoluchowski-like and then non-Smoluchowski regime 
sets in. If, however, t c <C t <C t q , then the evolution will be according to the 
classical equation. Notice, as h becomes small and the system is classical, t q 
will become large, which explains how (47) will reduce to the Smoluchowski 
equation. 

In this paper, we have given a formal proof, which is important for any 
new equation. It is our belief that examples will help in understanding 
(47) more. We wish to emphasise that the derivation does not assume any- 
thing special about the initial density operator (like, e.g., the Kubo-Martin- 
Schwinger condition). This generality is very important to note. However, 
we have employed ad hoc coarse-graining which, hopefully, does not destroy 
the novelty. 
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APPENDIX - A 



To prove (24), observe that 
d 



dE 



d_ 
~dE 

d 



J2(m\6(E-H)\m)-J2(n\d £t H5(E 

m ^ n 

-^j J2(n\d et H5(E - H) 
J2(n\d et H5'(E-H)\n) 



n 



d 



d(et) ^ 

as 



^(n\6(E - H)\n) 



d(et) 



The last equality follows because (n\m) = 5 nm implies that 
' ° n|) 5(E - H)\n) + (n\5(E - H) 1 ° 



\d(et) 

= {n\5(E - H)\n) + (n\5(E - H)\h 
= [(h\n) + (n\n)}5(E - E n ) = 0. 



d{et) 
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